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ABSTRACT 

The prototype system for Deep Space Network (DSN) clock synchronization by 
VLBI has been denionstrated to operate successfully over intercontinental base- 
lines in a series of experiments between Deep Space Stations at Madrid, Spain, 
and Goldstone, California. As predicted by analysis and short baseline demon- 
stration, the system achieves reliable synchronization between 2Gm and 64 m 
antenna stations with 17 and 37 K nominal system temperatures using under one 
million bits of data from each station. Semi- real-time operation is feasible since 
this small amount of data can be transmitted to JPL and processed within minutes. 
The system resolution is 50 to 400 ns, depending on the amount of data processed 
and the source intensity. The accuracy is believed to be comparable to the reso- 
lution, although it could be independently confirmed to only aLaut 5jus using 
LORAN C. 


INTRODUCTION 

The prototype for a semi-real-time system for DSN clock synchronization by 
radio interferometry was first demonstrated on a short baseline in August 
1972.^’^ A series of three experiments has now been conducted between Madrid, 
Spain and Goldstone, California, to demonstrate that the system performs as 
expected on intercontinental baselines. The series of experiments had three 
primary classes of objectives, all of which were achieved: 

1. To confirm that the system achieves the predicted resolution with the 
predicted amount of data; this implies that there arc radio sources avail- 
able which act as strong enough point sources over the long baseline. 

2. To check the system accuracy by demonstrating consistent results at 
different times of day and with different radio sources and, by direct 
comparison with other clock synchronization systems, to within the 
accuracy of these systems. 
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3. To gain experience in operating VLBI systems so as to uncover potential 
problem areas and facilitate the design of a higher-accuracy, operational 
c'ock synchronization system. 

Using one 64 m and one- 26 m antenna, several sources were available which were 
strong enough to achieve resolutions of from 200 to 400 ns with about 1 million 
bits of data and 50 to 100 ns with several million bits. VVe believe that the system 
accuracy is consistent with these resolutions. However, the absolute accuracy 
could be confirmed only to the approximate 5-10 ps accurac}' of LORAN C. 


DESCRIPTION OF EXPERIMENTS 

The three experiments were conducted on April 12, April 30, and June 11, 1973, 
On the first of these days, two 26m antenna stations were used, DSS 12 at 
Ck)ldstone and DSS 62 at Madrid. The last two experiments used the 64 m antenna 
at DSS 14, Goldstone together with DSS 62. To maximize the use of station time 
and to obtain a rough check on the accuracy, the last two experiments were run 
simultaneously with another VLBI experiment (conducted by J. Faneslow of JPL). 
The purpose of this other experiment was to measure the platform parameters, 
the source positions, and UTl. The two indei>endent experiments agreed in their 
estimates of clock offset to within about 10 ps, the approximate accuracy of the 
platform parameter experiment. 

The data were acquired and processed in the same manner as for the first clock 
synchronization demonstration. 3 Coherent bursts of about 0.31 Mbits of 
usable data, as limited by the computer memory size, were taken at 10 second 
intervals and recorded onto magnetic tape. Since the sampling rate was 
0.5 Mbps, the burst length was only about 0.62 second, but the burst frequency 
was limited by the 200 bpi magnetic tape density. Each full tape of data consisted 
of 72 bursts. However, some tapes contained less data due to tape failures. In 
the joint experiments, the source schedule of the platform parameter experiment 
was followed, which sometimes allowed only 2 or 3 minutes between sources. 
When only one tape unit was operational, the number of bursts of data on each 
tape was reduced due to the time required to change Lopes. Other short tapes 
were caused by tape errors. 

The tapes were mailed to JPL for processing. The maximum likelihood estima- 
tor function^ was calculated for each burst of data, and then the function values 
were accumulated in a nonoptimal manner over enough successive bursts of data 
to achieve an adequate estimator signal-to-noise ratio. The accumulation was 
nonoptimal only because fringe phase coherence was not maintained between 
bursts. This resulted in little degradation for strong input signal-to-noise cases. 
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when only 2 or 4 bursts were required, but in significant degradation for weak 
sources and when both stations had 26 m antennas. The system was not designed 
to operate in these situations because too much data is required for satisfactory 
operation. 


RESULTS 

Estimates of clock offset were obtained for two tapes of data on April 12, using 
two 26m antennas, and for ten and five tapes of data on April 30 and June 11, 
respectively, using one 26 m and one 64 m antenna. The nominal system tem- 
peratures were 17 K at the 64 m station and 37 K at the 26 m stations. Successful 
estimates were not obtained for all pairs of tapes on any day. On the first day, 
this was because the signal-to-noise ratios were too low on all but the strongest 
source. On the two days using the 64 m antenna, the failures are believed to be 
due mainly to tape failures. This is discussed further in the next section. 

The clock offset estimates are presented in Table 1, together with estimates of 
the resolution obtained and of the amoimt of data required for reliable measure- 
ment. The first three columns identify the case by date, time (GMT), and radio 
source name. The fourth column indicates the number of subcases into which 
the data were divided, and the fifth column indicates the number of coherent 
bursts of 0.31 Mbits of data that were combined noncoherently in each subcase. 
The sixth column is the average estimate of clock offset for the case, in micro- 
seconds, with positive values indicating that the clock at DSS 62 is early. The 
seventh and eighth columns are the estimated resolutions for one subcase and 
for the entire case in nanoseconds. For cases with more than one subcase, 
the resolution for one subcase is taken as the sample standard deviation, and the 
overall resolution is this standard deviation divided by the square root of the 
number of subcases. For the four caseshaving such alow signal-to-noise ratio 
that there was only one subcase, the estimated resolutions were obtained from 
the value of the estimator fimction. These estimates are quite unre’iable. The 
final column m Table 1 is the estimated amount of data required to obtain reliable 
estimates for that source, on that day, and with the system parameters used. 

It is the amount of data required for an estimator signal-to-noise ratio of 10 and 
a resulting resolution of 366ns. 

Data Requirements for Reliable Results 


Of primary importance is that reliable results were obtained with four radio 
sources with one million or fewer bits of data when using the 64 m antenna. These 
sources are 4C39.25, OJ287, DA193, and NRAQ190. The flux densities were 
not estimated directly from the data, because this requires considerable auJi- 
tiona! i^omputer time. However, from the estimator signal-to-noise ratios and 
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Table 1 


Clock Offset Estimates, Estimated -Resolution, n ' )ata Required 
for Reliable Performance 


Dat«. 1973 

Timf, GMT 

Source 

Vumber of 
subcases 

Number of 
bursts iRT 
subcase 

Kitunatts! 
CI«H k 4»lls4, t. 

Resolutions 
h>r one 
iubi'.ise, ns 

Resoliit-en 
for eetir • 
ns 

Minimum 
data re.,uired. 
Mbits 

4/12 

2132 

4C3925 

4 

36 

134 63 

345 

172 

(2 tapes) 

10 

4/30 

1756 


15 

4 

125.77 

340 

90 

1.1 


1812 

3C84 

8 

8 

125.17 

510 

180 

5 


1345 

D\V0224^67 

8 

8 

124 41 

470 

170 

4 


1906 

NRAyiOO 

6 

8 

125 59 

220 

90 

09 


2118 

3C 171 

i 

38 

1227 

400* 

400- 

14 


221'3 

4C39.25 

H 

4 

12549 

190 

70 

04 


2242 

Q1JS53 

1 

72 

126.6 

3<‘0* 

370« 

22 


23(K> 

Oj:8T 

23 

2 

125.56 

240 

50 

0.3 

5/1 

0007 

3C37I 

I 

SI 

123 3 

221> 

220- 

6 


0022 

4C3925 

9 

4 

125 24 

150 

50 

0.2 

6/11 

1726 

IM193 

16 

4 

2ai 16 

220 

60 

0.5 


1S48 

OI3G3 

S 

8 

203.39 

470 

17u 

4 


1905 

DA 193 

16 

4 

202 71 

190 

50 

04 


1923 

4f:.39.2.5 

16 

4 

203 25 

210 

50 

04 


1911 

3r;ri 

I 

71 

2(M» S 

1000 ■* 

1000* 

50 

*Kor these 

cas€^ the si;^ri.i!-tj 

)-noi'.4* ratio is so 

low th.it the 

esbmates of resolution am! data 

rtv^uirtal m.iV h; 

»vr I.oce er.ors. 


the nominal system temperatures, vve estimate the correlated flux densities for 
these sources to raf^ge from 1 to 2 f. u. There are enough radio sources of this 
strength that clock synchronization is attainable at any time of day between DSS 
14 and Spain or Australia. Similar performance will be achie' able between Spain 
and South Africa when DSS 63, a 64 m station, is operational in Spain. 

Resolution and Consistency 


The results for the ten cases of April 30 are shown as a unction of timeinPigui'e 1. 
The straight line has a slope equal to the rate of change in the clocks,— 

0. 6ps/day, as predicted by the fringe rate data b.ken in the platform parameter 
VLBI experiment on that day. The error bars shown are the estimated resolu- 
tions. At best, the accuracy would be the resolution plus a degradation of up to 
about 100ns due to propagation and geometrical effects. This should be added 
in a mean square sense. (There is also a constant error due to time delays in 
the stations, which would not show up in the results.) 

Overall, all of the clock offset estimates fall as close to the straight line as ex- 
pected except for the case for DW0224+67 and the two cases for 3C371. Al- 
though these sources are quite weak, especially 3C371, there appears to be 
some significanterror for these cases. The cause for this error has notyet been 
satisfactorily explained. However, one likely cause is that the tim.e delay pre- 
dictions were in error due to the high declinations of the sources, +67deg. for 
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Figure 1. Clock Ctffset Estimates for April 30, 1973, with lo Resolutions 

D\V0224-!^«17 and +69deg, for CC371. The highest declinations for the other 
sou»-ces are -oOdcg. forOL.'>;73 and -41 deg. for 3C84. Furthermore, the geom- 
etry was so unusual for 3C371 that the received frerjuency was high ■'r for DSS C2 
than for DSS 14, just the opposite from the usual case, be< ause the antennas were 
pointed on the opposite from usual sides of due ; orth. This can happen only with 
high declination sources. 

Comparison to LORAN C and Platform E?q)eriment 

The clock offset estimates obtained -n these experiments agree with estimates 
made on each day using LCRAN C, to within the approximate 5-lOps accuracy 
lisat could be expected. For the last two days, the results also agree with the 
v'LBI platform parameter experiment, totheaccuracy of that experiment. There 
was, however, a constant 200ps offset with respect to LORAN C. This was 
determined to arise because some of the 1 second period reference signals at 
D5>3 <52 occurred 200 ps early oecause of reference to the wrong edge of a 200 ps 
duratior mjlse. 

Tne comparisons to LORAN C and the platform experiment are shown in Table 2. 
The first column lae ifies the experiment date. The second column is the 
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Table 2 


Comparison of Clock Offset Estimates 


iw. 

I9T3 


U»K\\ <• I.OKW r. 

|»NStC* i:\irtr\tp >r,\ Oss 14 
i . MS tol.t»RA\<; 


LORAN i 

*^tanjtr 1 >^S ♦li t4» 
IKS 12 U MS 


jt^ramrtrr 

•> 


< hK r«liRMt ' of 

llrf\ fipmttM M. 


4 12 62 -1 -61 - 65 

4 ?0 -74 *7 67 64 -74 

6 11 *1 4 - I 9 *3 


difference between the DSS 6? and the LORAN C Mediterranean chain clock. The 
nominal accuracy is about 2iis, but ambiguities oi.en cause errors in 10/is incre- 
ments. TTiif third column is the sum of all knovm corrections from the LORAN C 
Mediterranean clock to the DSS 12 or 14 clock, that is LORAN C to the U.S. 
Naval Observatory (USND), USNO to the National Bureau of Standards (NBS), 
NBS to the Gcldstone master clock and the Goldstone clock to DSS 12 or 14. Each 
of these legs introduces some error, typically on the order of 1 ps. The fourth 
column is the estimate of clock offset between the stations obtained from LORA N C . 
The last twc- columns are the offsets obtained by the VLBI platform parameter 
ejqjeriment and by this experiment, corrected for the known 200 ps error at DSS 62. 

The clock offsets measured in tiiis experiment agree with LORAN C to within an 
average of 3ps and to the VLBI platform experiment to an average of lips. 
These results are consistent with the accuracies of LORAN C and the platform 
experiment. 


PROBLEMS, CAUSES, AND SOLtHTIONS 

The two major problems in the series of experiments were that results were not 
obtained for any day until the data from the last experiment were processed and 
that good results were not obtained at all for a number of cases for which the 
signal-to-noise ratio should have been adequate. A number of possible causes 
for lack of results were studied while an unsuccessful attempt w~s being made 
to process the data. These investigations will influence the design of the final 
operational system. 

The ac'ay in achieving results was caused primarily by the 200ps error in the 
predicted clock offset. This area of clock offset was not searched until after 
the large offset error was discovered by the VLBI platform parameter experi- 
ment ''rs. Even then, there wa^ . delay in obtaining good results due to a dis- 
crepancy between the s^'fware and the hardware which was caused by a system 
modificat'on made after the short baseline experiments. One of the ramifications 
of this ^ was overlooked. Aside from leading to a more complete analysis 
of die system, this processing problem will have no effect on the future system. 
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Search Range and Variable Bandwidths 

The 200ns offset error will have an effect on the future system. It pointedly 
illustrates that large errors can occur, even when the system is siq^sedly cal- 
ibrated with LORAN C. It also illustrates the problems which could occur if 
synchronization is lost and the best available timing information available is by 
radio bn^adcast, with errors on the order of milliseconds. An opeiational sys- 
tem should be capable of efficiently searching as wide a region of clock offsets 
as possible. This is best accomplished 1^ enabling data acquisition with a rum- 
ber of selectable bandwidths and sampling rates. Lower bandwidths and sam- 
pling rates enable wider searches because the wide region encompasses fewer 
samples. Programmable digital filtering would be the most flexible and econom- 
ical method for achieving this. 

Local Oscillator Frequency Adjustment 

One possible cause for bad data w’as that the third (10 MHz) local oscillator (LO) 
was offset at the Goldstone station in order to compensate for t>'e fringe rate or 
doppler frequency differ**- 'e. This LO was manually adjusted for each case, 
and it is possible that the frequency synthesizer was improperly set for some 
cases. This problem can be avoided by accounting for the fringe rate «itirely 
in the data processing. There would be &ome cost in computer time. 

S 3 mchronization of Samples to Station Clock 


Anoffier possible trouble area is in s}mchronization of the samples to the station 
clocks. these experiments, the hardwc«e was resynchronized to the frequency 
and timing system (FTS) Ipps and 1 MHz signals at the beginning of each burst 
of data, i.e. , every 10 seconds. It would probably be more reliable to sjmchro- 
nize the system to the station clock only once for any experiment. Since the reso- 
lution of the present FTS is only on the order of 0. 1 to l.O^iS, a higher- resolution 
clock will be required in any case for the more accurate clock synchronization 
system. This high- resolution clock could either be a part of the clock S 5 mchro- 
nization hardware or be in the form of a new timing system. 


CONCLUSION 

The prototype VI.BI clock synchronization system has been demonstrated to 
operate as predicted over transcontinental baselines. Resolutions of 50-400 ns 
were achieved, and the accuracy is believed to be approximately the same. A 
wideband system is being designed which will achieve a resolution of 1-10 ns, 
depending on the available receiver bandwidth. The accuracy of this system 
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will be limited by atmospheric and geometric effects and bj' errors in calibrating 
time delays in the stations rather than by the system bandwidth and resolution. 
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QI LSTIOX AND ANSWKR PKRIOD 


MR. CHf: 

Are there any questions? 

DR. WINKLER: 

What clocks have you used between the two stations for which you showed us 
before the last slide, the second to the last slide? 

DR. HURD: 

Well, all the experiments used rubidium. I think there might have been a maser 
at one station, but there was at least one rubidium in the system. That is the 
advantage to this tjpe of system; since we are taking the data over such a short 
period of time, the oscillators ability isn't a problem. 

OR. WINKLER: 

May we see the second to the last slide again, Ix'cause I have several questions 
on that. 

(Slide, i 

DR. AINUI.ER: 

There is something on that slide which, I think, even ’f you have two very f>oor 
rubidium standards is simply incomprehensible. 

DP.. HURD: 

Well, to fill this buffer in the existing svscem takes six tenths of a second — it 
is a 320 kilobit buffer being filled a': a 500 kilobit rate. 

DR. WINKLER: 

Yes. 

DR. HLHD: 


So, to — 
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DR. WINKLER: 


But my question aims at clarification of the meaning of the error bars which you 
have here which evidently must be connected with your estimate of internal 
precision of each individual synchronization experiment. Evidently there must 
be very much larger external errors. 

There simply is no clock which v\tmld jump around four, five or six micro- 
seconds between one and two hour intervals. Even the worst crystal would do 
much better than that. 

So, you must have considerable external errors, which are completely ignored 
here. 

DR. HUKxJ: 


Well, as I tried to point out, ^ haven't explained why the few larger error residts 
ocr- f-red. These were weaker sources; they were in peculiar positions in t*- 
sky; and we didn't calibrate the atmosphere. The atmosphere couldn't cause 
that much error, though. 

DR. WINKLER: 

What sense do the error bars make then? If you put error bars on which appear 
to have absoluteh' no significance whatsoever — 

DR. Him.- 

Oh, the error bars represent the errors due to statistical errors and due to 
receiver noise. 

The error bars are the ctjseived standard deviations of the errors for that 
source. In other words, they are how consistent the measurements were for 
that particular source. For example, for the first point on the curve, I believe 
that one tape of data was divided into 18 cases. Each case was an accumulation 
of four of these batches of data. So that we had for aach batch of data roughly 
one million bits of data for each of 18 cases. ITie error bars represent the 
standard deviation of the estimate for those 18 cases. In other words, it is how 
well we would have done, you know, if the system had been perfectly calibrated. 

MR. CHI: 

The error bars enumerate the number of data points ? 
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I)K. BI RD: 


Well, no. If you gather statistical data and you want to make an estimate of it 
as to how accui'ate you are estimating some parameters, you divide it up into 
several batches, and compute the sigma, the standard deviation of these 
batches, and then you divide that sigma by the square root of the total number 
of batches to get the standard deviation on the entire amount of daU, and that 
represents the error just due to receiver noise, not due to atmospherics, not 
due to anomalies in the clock changing during the time. Well, you can see, 
even with the worst cases, there are not verj- man3’’ signals from the line, four 
or five sigma. 

DR. WINKLER; 

Yes, but what you seem to claim is that your clock essentialh- has shifted by 
several micrr>seconds every few hours. 

DR. HERD: 

No, I am claiming that that measurement was a bad measurement. I am claim- 
ing that it was a poor measurement, it is a bad result. I am not just showing 
the good results. 

MR. CIII: 

There is a question. Dr. Reder. 

DR. REDER; 

Isn't there realh’ a misunderstanding between the two of you? Are these 
measurements taken on the same da\* ? 

DR. HURD: 


Yes. 


DR. REDER: 

They are taken on the same day? 
DR. Wl vTCLER: 

Yes. 
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DR. REDER: 


Then I vvill be included in the misunderstanding. 


DR. HITID: 

These are results on the same day, and they don't all agree. The results for 
the stronger sources do agree, as well as we expected. The results from some 
weaker sources have unexplained large errors. 

DR, ClARK: 

Tom Clark, NASA-Goddard. I have tw-o questions or comments. 

I would first of all be curious in making these clock calculations which you had 
here. You obvioush' have made some assumptions as to the positions of the 
sources, and the length of the baseline. That is, those factors included in the 
geometry. 

It would be interesting to know what assumptions were made in that case. 

Second of all, for the much vvider band width system which you described, there 
is one effect which you didn't indicate how it would be taken out, about which 1 
am curious. This is differential phase response over the pass band, which is 
going to cause an instrumental delay dispersion. 

Do you plan some sort of instnunental phase calibration to be taken v.ith this 
data? 

DR. HURD: 

I feel that since I am using the full continuous band width, phase stability' 
requirements on the system will be less critical. 

On the second point of view, to answer your question as to how good I expeet 
this to be; the projections are for a total system stability, from and including 
the front end maser and the receiver system, which is supposed to be stable to 
Id centimeters, over a 12 hour period or so, with ceiiain assumptions on how- 
much the temperature can change and things like that. 

So, at the 10 nanosecond level, I think that tne system stabilitj- will not play an 
important role here. 1 think that atmospherics -.vill be dominating. 1 think that 
how good you do with the system depends on wh.it you do with the atmosphere. 
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To get 10 nanoseconds, I think you would agree that we can probably a priori 
noodel the atmosphere that well, especially at S band. 


DR. CLARK: 

The question I was asking, however, had to do with the differential phase 
response over the pass band. You are going to just have to rely on the networks 
themselves calibrating that out for you. 

DR. HURD: 

The entire differential group delay is 10 centimeters, it can vai-y 10 centimeter's 
over a day. 

DR. CLARK: 

No, I am talking about at one instant of time, over the pass band, the group 
delay wliich you are measuring, the time delay number you are measuring is 
essentially rate of phase change over the pass band. Now, from that, that number 
has to be subtracted from the instrumental phase response, and I was asking — 

DR. HURD: 

1 don't even measure the phase. 

^R. CLAK... 

Well, that still has to be defined. This is a group delay, and group delay is 
A0/Acj , so the instrarmental phase response, the dispersion of the filters in the 
pass band are an important factor to you. 

DR. HURD: 

The tentative spec on the system is 10 centimeters for eight hours total change. 
That is, change in the total delay through the system. 

Clearly, the differential delay at the two edges of the pass band is going to be 
much better than the total delay in the whole r , stem. 

DR. FLIEGEL; 

Henry Fliegel, JPL. 
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1 believe the question w-as that if you want to measure the absolute difference 
between two clocks at tw'o different stations, that you have to calibrate the 
different phase delays and different channels that you are using. 


DR. HURD: 

I am not using different channels. 

DR. FLIEGEL: 

Or across the one channel that you are using then, if it is broad banded. 

DR. HURD: 

Bj' utilizing the entire continuous oand width, what you are buying, in a loose 
sense, is the average change in the group delay as a function of frequency, 
rather than being particularly concerned with the differential ,'hase delay at two 
edge bands, and I think it is clearly easier to average this thing out. 

DR. FLIEGEL: 

Let's get together over coffee. I think that there is an issue here that we are not 
getting across. 

DR. HURD: 

That is not my area of expertise, in any case. 


MR. KAUFMANN: 


I have a question in regard to this system relative to moon bounce. Is it more 
economical, what kind of savings do you envision, and is the performance 
greater, and if so, how much better? 

DR. HURD: 

Well, performance is cetainly much greater. The moon bounce system 
accuracy is 10 to 20 microseconds. 

The major cost of the moon bounce system is in maintaining the system. There 
is over a million dollars worth of equipment involved in the system, and if you 
make reasonable estimates as to what your average long term maintenance 
figures are going to be for this amour.t of equipment, it turns out to be quite a 
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lot of money. The cost is considerably more than the amount of money which is 
now budgeted for maintaining that system, which is 40 thousand a year. 


1 mentioned to Bill that a VLBI system would cost in the neighborhood of a 
counle hundred thousand to build. This includes three data acquisition systems 
for a reliable digital system. 

MR. CHI: 

Thank you. Bill. 
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